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Abstract 

A previously tested differential equation method for generating low temperature series 
expansion for diagonal spin-spin correlation functions in the d = 2 Ising model is extended 
to generate the non-universal terms for arbitrary separation of the spins. This extends 
the earlier calculations of these correlation functions. 

The two-dimensional Ising model is one of the very few exactly solu- 
ble systems Many thermodynamic functions for this system have been 
calculated (exactly or approximately) for this system in classic papers 
long time ago^^'^l Among the interesting objects of this model are the 
(static) spin-spin correlation functions < doo'^rnn > W't^^'t^lA few years 
ago, the author in a collaboration , used a differential equation approach 
to the diagonal correlation functions Sn =< o-Qoann > of this model ob- 
taining the solutions as high and low temperature series expansions 
in suitable expansion variables^^l Although the high temperature ex- 
pansions generated this way were very successful in getting the exact 
coefficients for arbitrary n,the low temperature expansion was not very 
successful as it was realised that these expansions have a universal part 
which comes from the series for the exact expression for the sponta- 
neous magnetisation of the system, and the rest of the expansions have 
a non-universal character i.e. they depend on the particular value of 
n.In that paper we did not say anything about how to go about cal- 
culating these non-universal terms. Soon afterwards, the same authors 
calculated the exact expressions for these correlation functions for val- 
ues of n upto n = 6 However these expressions become more and 
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more complicated as the value of n grows and because of this , these 
were given only for values of n upto n = 5 in the publication at that 
time. It is clear that if we want to write down the expressions exactly 
for large values of n the computation time and the labour involved 
starts becoming unmanageable. Hence it seems worthwhile to look for 
a method for calculating these non-universal terms in the low temper- 
ature series expansion. At the same time ,to the author's knowledge, 
there exists no systematic method for calculating the low temperature 
expansion for these correlation functions. The purpose of this paper is 
to show how we can calculate these very terms by the same differential 
equation approach that we had presented earlier. For the exact coef- 
ficients for the high temperature series for these correlation functions 
the reader is referred to[7]. In this short communication we follow the 
conventions adopted in our earlier paper to avoid any confusion in the 
results. 

We take the Hamiltonian for the model to be 

H = - iJm,j(^i+i,j + J2(Ti,j(^i,j+i) (1) 

where Ji and J2 are the exchange constants of the model. Let/3 = 
{UbT)'^ and as usual define 

k> = sinh{2(3Ji)sinh{2(3J2) (2) 



applicable for T > Tc, and 

k< = k;' (3) 

applicable for T < where Tc is the critical temperature of the 
model, which is given by 

A:>(/3c) = k4(3c) = 1 (4) 
The analysis then proceeds as follows. The variable t is defined as 

=p are defined to be < (jQQCJnn > for T < and T > Tc respectively 
and we define two related functions by 

a^^_ = tit-l)j^lnSn,--\, (6) 



d 1 

an,+ = t{t - l)—lnSn,+ - -t. (7) 

Jimbo and Miwa showed that (7„ = (Jn,± (t) satisfy the ordinary differ- 
ential equation!^] 



[t{t - lyX-n' lit - 1)< - it - 1)< - - 4 



(8) 

It was shown in ref . [7] that if we consider the high temperature phase 
then we can generate the series expansion for the spin-spin correlation 
functions rather easily, as there are graphical methods available for these 
expansions which give us the first term in the series with a little bit 
of combinatorics and then many terms of the series can be generated 
recursively using the differential equation. The program was very suc- 
cessful when carried out this way. However when the same equation was 
applied to the low temperature phase of the theory for a series expan- 
sion,beginning with the first term as 1 (which is the appropriate term 
for zero temperature) it ran into a problem. The recursive procedure 
gave us an expansion 

= i: (-!)'( I )fcJ'+o(fci'+^) (9) 

where ((^)is the binomial coefficient in the expansion of the square of 
spontaneous magnetisation, given by 

= lim < (7QQ(7mn >= (1 " k^)'^ (10) 

and the terms of the order k'^'^'^ in equation (9) depend on the specific 
value of n. Because of some reasons , the calculation of the n dependent 
terms could not be carried out at that time. Although we think that 
they carry important information about the detailed nature of spin-spin 
correlations along the diagonal. 

The reason for our failure at that time was the fact that once we 
take the series and substitute it into the equation the coefficients of 
the terms of the order k<^^^^ and /c<''^^^ get coupled together in the 
same term along with the other coefficients. This makes it impossible to 
determine the first n dependent term in the series without recourse to 
some other information. 



Fortunately a closer examination after all these years has revealed 
that the first n dependent term in the low temperature expansion of 
Sn,-is given by the expression 

Where Sn- is given by 

oo 

This relation ^ has been tested for n < 6 and is conjectured to hold 
for arbitrary n. Armed with this information we go on to calculate the 
series for Sn - for values of n upto 10. It is seen that with this value of 
the coefficient, the series can be generated term by term, and this can 
not be done for the other possible values of the coefficient that we can 
try.However the most important criterion for choosing this function as 
the correct value of the above term is the fact that it gives the correct 
coefficients for all the subsequent terms of the series for known exact 
correlation functions. We carry out the expansion for the first thirty 
terms in series from n = 6ton = 10. The expressions for these 
functions upto n = 5 have been omitted because the exact expressions 
for these are rather simple in this case, and are already available in the 
literature. 

Cn,m = 0, for m = odd no. (13) 

n = 6 



^This relation seems to imply the equality (^p-)"^^ ('^'i,- + <S'ra,-)|fe=o = (gp')""'^^ M^|fe=o-J-H.H. 
Perk has informed mc that the proof of this equality is easy using the quadratic recurrence relations 
for the correlation functions given in the ref.[10]. 



C6,14 — 



62755 
4194304 



Ce.ie — 



5885649 
536870912 



n = 7 



n = 8 



n = 9 



C6,18 = 
C6,22 = 
C6,26 = 

Ce.so = 

^7,16 '- 
^7,20 '- 
^7,24 = 
^7,28 - 

C8,18 = 
C8,22 = 
C8,26 = 
C8,30 = 

C9,20 
C9,24 
^9,28 



18102051 

■ 2147483648 

1502308899 
274877906944 

135253617001 

■ 35184372088832 

12860686137607 

■ 4503599627370496 



13852377 



1073741824 



529859407 
" 68719476736 



91202314827 
17592186044416 



1052329791057 
281474976710656 



48531703 



4294967296 



973272573 
137438953472 



172699772293 
■ 35184372088832 



4071839681233 
1125899906842624 



688829933 
68719476736 



114906443431 
17592186044416 



1307717138979 
281474976710656 



C6,20 = 
C6,24 = 
C6,28 = 



^7,18 = 
C7,22 = 
^7,26 = 
^7,30 = 

C8,20 = 
C8,24 = 
C8,28 = 



^9,22 = 



^9,26 — 



^9,30 — 



115143369 
■ 17179869184 



39992200521 
" 8796093022208 



927283251903 
281474976710656 



41988023 
" 4294967296 



1719298113 
■ 274877906944 



307615754733 
70368744177664 



1821460877687 
562949953421312 



302177557 
■ 34359738368 



25689376531 
"4398046511104 



1178557994739 
■ 281474976710656 



2191660887 
274877906944 



384494320401 
70368744177664 



281692820385 
70368744177664 



n = 10 



Cl0,22 — 



2470054587 
274877906944 



Cl0,24 — 



64020065109 
8796093022208 



Cl0,26 — 



212856341061 
■ 35184372088832 



Cl0,28 — 



1441790227059 
■ 281474976710656 



Cl0,30 = 



19816146315693 



4503599627370496 

We have stopped at the terms of the order in all these expansions. 
But if needed these series can be continued to higher powers of A;<. 

Thus in this paper we have given a scheme of calculating the low 
temperature expansions of the diagonal correlation functions in the 
two-dimensional Ising model and we hope that this method will be 
useful in getting information about similar quantities in other models 
for which one can write down differential equations. 

I am thankful to J.H.H.Perk for his comments and suggestions on 
the first version of this note. I am also thankful to him for bringing to 
my notice several recent references including ref. [5],[6] and [10]. The 
ref. [5] is specially important, for the series for susceptibility has been 
calculated in there. The incompleteness of the list of references is, 
however, due to the author's unfamiliarity with the later works. 
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